Thermopower of Aharonov-Bohm Interferometer: Theoretical Studies of Quantum 

Dots in the Kondo Regime 
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We report on the thermopower of an Aharonov-Bohm interferometer (AB) with a quantum dot 
in the Kondo regime. The thermopower is anomalously enhanced due to the Kondo effect as in 
heavy fermion systems. In contrast to bulk systems, the sign of the thermopower can be changed 
by adjusting the energy level scheme or the particle-hole asymmetry of a dot with the gate voltage. 
Further the magnitude and even the sign of the thermopower in the AB ring can be changed at will 
with varying either magnetic fields or the gate voltages. 
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Observations!^, |2| ^ of the Kondo effect in quan- 
tum dots have opened a testing ground of the quantum 
effects of electron wave functions and many-body effects. 
In contrast to bulk systems where the phase is washed out 
by several scattering processes, the phase coherence of an 
electron's wave functions can be preserved in nanoscopic 
systems. A typical system to measure the phase co- 
herence is the Aharonov-Bohm(AB) interferometer. Re- 
cently a phase-sensitive Fano resonance was observed jsj 
in the differential conductance for a quantum dot system 
which contains both a resonant current path and a direct 
one. 

An asymmetrical differential conductance suggests 
that the thermopower may be the right experimental 
tool for the study of the phase-sensitive Aharonov-Bohm 
interferometer. There are indeed some studies of ther- 
mopower in AB ring geometries |^ for noninteract- 
ing electrons. The thermopower of electrons is sensi- 
tive to the particle-hole asymmetry in the density of 
states(DOS) and the energy dependence of the electron 
scattering rate. In heavy fermion systems(HFS), the 
enhanced thermopower and its sign are determined by 
the energy dependence of the Kondo resonance scatter- 
ing rate of the conduction electrons off the magnetic 
ions^, H]. In mesoscopic systems, the thermopower will 
probe sensitively the asymmetrical shape of the trans- 
mission probability. 

In this paper, we study theoretically the thermopower 
of quantum dot systems in the Kondo regime. The ther- 
mopower of a quantum dot is anomalously enhanced due 
to the Kondo effect as in the HFS. In contrast to the 
HFS, the energy level scheme or the particle-hole asym- 
metry of a dot can be modulated continuously by the gate 
voltage capacitatively coupled to the dot and the sign of 
the thermopower can be changed from negative to posi- 
tive. When a quantum dot is inserted in the AB ring(see 
Fig. the Fano interference leads to a more dra- 
matic change in thermopower. Due to the Kondo effect 
of a quantum dot, a new Kondo-resonant current path 
opens below the Kondo temperature and interferes with 



the direct tunneling path leading to the Fano interference 
in the transmission probability. Adjusting both the AB 
phase by varying the magnetic field and the tunneling 
matrices by changing the gate voltages, the transmission 
spectral function can be controlled to take several differ- 
ent shapes. In addition to the differential conductance, 
the thermopower can probe the shape of the transmis- 
sion probability since the thermopower is sensitive to the 
particle-hole asymmetry in the transmission probability. 
Combining the Kondo effect and the Fano interference 
leads to an enhanced thermopower of the order of fcs/e. 
The magnitude and even the sign of the thermopower 
can be controlled with varying the AB phase. 



We may describe the quantum dot using the Anderson 
impurity model, Hd = '^dJ2a=t,i'^i'^a + Un^n^, when 
the number of electrons (A'^) in a quantum dot is odd and 
the spin of the highest-lying electron is unpaired. Here 
6(1 is the energy level of the highest-lying electron with 
unpaired spin in a quantum dot and U is the Coulomb 
interaction. The left and right leads are described by the 
noninteracting Hamiltonian H„ = y^,r e rc^- c - with 

o P ^ka pk pka pka 

the lead index p — L,R. The tunneling of electrons be- 
tween two leads via the direct tunneling(Tifl) or through 
a quantum dot(Vdp = V*d) is described by the Hamil- 

tonian, H, = [imEkk' [Tlr c[^,^c^-^ -\- H.c] + 

i^/W] Eka [^piikJc + H.c] . The AB phase 

(j) — 27r$/<I'o is contained implicitly in the tunneling ma- 
trices in a manner that VdhTLRVRd = iVdiTtfiVRdle**. 
$ is the magnetic flux passing through the system as 
shown in Fig. |I] and $o — hc/e is the flux quantum. 
Recently, this model system with one conduction chan- 
nel was studied for the conductance using the equation 
of motion method |ll| and the numerical renormalization 
group in equilibrium! 1 4. 



The electric and heat current operators can be de- 
fined as a change in the number of electrons and the 
total energy per unit time in the left electrode, II = 
e[NL,H]/ih and Ql = —[Hl, H]/ih, respectively. Here 
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Nl = yir - c - is the number operator and is 

^-~'t^OL Lka Lka 

the Hamiltonian of the left lead. We use the Keldysh 
Green's function method p3|, Q to write the electric and 
heat currents. Summing over all the multiple tunnelings 
between two leads and using current conservation in a 
steady state, the electric and heat currents can be ex- 
pressed in terms of the Green's function of a dotp, R|. 
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The numerical factor 2 accounts for the two spin direc- 
tions and fp{w) — fito — fJ-p) is the Fermi-Dirac thermal 
distribution function of p = L, R lead with /ip = —eVp. 
The transmission spectral function T{lo) is given by the 
following expression in the wide conduction band limit 



T{uj) = To + 2r^gn{l ~ To) cos^ReGJ 

+r [To -<?(!- To cos^c/.)] ImG^. (2) 



Here is the retarded Green's function of a quantum 
dot and F = (Fl+F7j)/(1+7) is the Anderson hybridiza- 
tion. Fp — TrNp\Vdp\'^ measures the hopping rate of elec- 
trons between the quantum dot and the leads, where Np 
is the density of states of lead p = L,R. Tq = 47/(1 + 7)^ 
is the transmission probability due to the direct tunnel- 
ing, and 7 = tt'^Ni^NjjITi^jiI'^ is the dimensionless mea- 
sure of direct tunneling of electrons between the two 
leads, g = 4F/^Ffl/(Fi -I- Tji)^ is the maximum dimen- 
sionless linear conductance through a quantum dot in the 
absence of the direct tunneling and also measures asym- 
metry in the coupling of a quantum dot to the left and 
right reservoirs. The interference effect is included in the 
second and third terms of T{lo). The above expression 
for T{uj) agrees with Eq. (2) in Rcf. 

We compute the dot's Green's function Gd using the 
noncrossing approximation(NCA) for the infinite U An- 
derson model in the Kondo regime. The NCA||l^, ^ 
p^ , [ig|] has been successfully used for the study of the An- 
derson model except for the nonanalytic behavior pnlat 
a temperature far below the Kondo temperature TR-. plfl 
In the NCA self-energy equations, the renormalized An- 
derson hybridization should be used instead of the bare 
Anderson hybridization p2| . The multiple tunnelings be- 
tween the two leads result in the flux-dependent renor- 
malized Anderson hybridizations which in a wide con- 
duction band limit are given by the equations. 
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(1+7)2 



'iisin0 (3) 



In equilibrium, the two thermal functions are equivalent 
or = fn. so that the total Anderson hybridization 
F = Fi + Tb, is independent of the AB phase and 
so is GJJ. T{uj) also remains invariant under the inver- 
sion of the magnetic flux: $ ^ — $. The Onsager re- 
lation T{uj, —(j)) = T(a;, 4>) in equilibrium becomes bro- 
ken when a finite source-drain bias voltage is applied 



(/i Ir)- The Kondo temperature Tr- is indepen- 
dent of <j) and can be estimated in the J7 — > 00 limit 
by the equation Tk = Dy'N{0)J exp (-l/iV(0) J) with 
N{0)J — 2F/7r|erf|. Note that the direct tunneling sup- 
presses the Kondo effect. 

At high temperature above Tk, the current flow 
through a quantum dot is blocked due to the strong 
Coulomb repulsion (Coulomb blockade). Electrons flow 
from the left reservoir to the right one only via the di- 
rect tunneling. With decreasing temperature below Tk, 
the Kondo resonance peak at the quantum dot develops 
close to the Fermi level. The newly opened current path 
interferes with the direct tunneling path. This Fano in- 
terference transforms T(a;) into various different shapes 
depending on the AB phase (j). The general structure 
of T{uj) near cu — can be read off from the Eq. (||). 
It is well known that — ImGJ develops the Kondo reso- 
nance peak with its width of the order of Tk near uj ~ 
while ReG^ varies very rapidly over the energy scale of 
Tk near lu — with a dip just below to = and a peak 
above oj ~ 0|^. The overall shape of the transmission 
spectral function T^u) is determined by the value of the 
AB phase and the sign of Ac [see Eq.(0)], 



Ac = ro-ff(l-Tocos2 



(4) 



A typical Fano interference pattern consisting of a dip 
and peak structure is expected when cos</) 7^ 0. At 
coscf) — 0, T{uj) has a dip(peak) resonance structure if 
Ac > (<)0, respectively. 

In our numerical NCA work we consider a symmetri- 
cally coupled dot {g — 1 or Tl = Tr) with the energy 
level scheme(ELS): the iV — 1 state lies lower in energy 
than the N + 1 state, En-i ^ En+i, where N is the 
number of electrons in a dot. Other ELS's will be dis- 
cussed later. In the Anderson model picture, E^^i = 
(empty). En = erf (singly occupied), and E^+i — 2ed+U 
(doubly occupied). In practice, we take the limit of 
U — > 00 1^ and the DOS of two leads is assumed to 
be Lorentzian of bandwidth D. 

For the weak direct tunneling, the Fano interference re- 
mains weak and the Kondo-related peak persists in T(uj) 
over all (f>^^- The AB phase 4> dependence of T{uj) be- 
comes stronger with increasing direct tunneling ampli- 
tude To- The case of Tq = 0.5 is displayed in Fig. || and 
the shape of T{uj) near w = is strongly sensitive to the 
AB phase 4>. Since Ac = at = 0° and 180°, T{uj) has 
the dip-and-peak structure of ReGJJ (see Fig. ||(a) and 
(c)). 

When coscj) = {(j) = 90°, 270°), the spectral shape of 
T(cj) is wholly determined by both ImG^ and the sign 
of Ac. Ac remains always negative for = 1 so that 
the Kondo resonance peak persists in T(uj) for the whole 
range of ro(0 <JTo < !)• The Fermi-hquid relation^, 
ImGJJ(O) = — 1/F at T = OK, leads to one exact relation: 
r(0) = g for any interacting dots. T(lo) reaches its maxi- 
mum possible value g = latuj = 0,a unitary Kondo res- 
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onance tunneling. When the background direct tunnehng 
amplitude Tq is increased from to 1, |Ac| approaches 
zero and the Kondo-related peak becomes smaller 

Except for close to cf) = 90°, T{lu) of asymmetrically 
coupled dots ((7 < 1) is not much different from the sym- 
metrically coupled dots|^. Since Ac = Tq — .g at = 90° 
can change its sign, the Kondo-related peak in T{u!) 
transforms into a dip at = as the value of Tq crosses 
g. T{lu) is flat near the Fermi level when Tq ~ gp2[. 

The thermopower of a quantum dot in a two-terminal 
configuration can be found in an open circuit (/ = 0) by 
measuring the induced voltage drop across a quantum 
dot when a temperature difference between the two leads 
is applied. The thermopower S is defined by the relation, 



S — — lim 



Vl - Vr 



(5) 



Expanding the expressions for / and Q in Eq.(^ up to 
the hnear terms of SV = Vl - Vr and ST ^ Tl - Tr, 
the transport coefficients can be expressed in terms of 
the integral, /„(T) = [2/h] J duj uj'^T{uj) [-dj/duj]. I = 
LiiSV + Li2(5r and Q = L21SV + L226T where Ln = 
e2/o(T), L21 = L12T = -e/i(r), and L22 = h{T)/T. 

The thermopower S = —Ii/eTIg probes the particle- 
hole asymmetrical part of T{uj). To begin we compute 
the thermopower in the absence of direct tunneling or 
To = 0. In our choice ofthe ELS in a dot, E^-i «C E^+i, 
the Kondo resonance peak has more spectral weight on 
the electron excitations. Since the electron excitations 
are the main carriers of charge and heat, the sign of the 
thermopower is negative. From the study of the heavy 
fcrmion systems (HFS), it is well known that the ther- 
mopower is anomalously enhanced due to the Kondo 
effect and is of the order of fcB/e(« 86.17/j,V/K) near 
T = Tk In normal metals, the thermopower is 

of order /iV/K. The thermopower is also enhanced in a 
quantum dot in the Kondo regime, as shown in Fig. ^by 
the solid line. 

To study the effect of the Fano interference on S, we 
now turn on the direct tunneling. The computed ther- 
mopower S{T) for To = 0.5 is displayed in Fig. ^ for dif- 
ferent AB phase (f). When </> = 0° , the electron transmis- 
sion is high while the hole transmission is low[Fig. ||(a)]. 
Since electron excitations are the main carriers of charge 
and heat, the thermopower is negative and its magni- 
tude is of the order of kB/e near T = Tk due to the 
Kondo effect. With increasing AB phase, the dip-peak 
structure in Tiuj) transforms into a Kondo-related peak 
at cj) = 90° [Fig. 1(b)] . In this case, T(w) is more or 
less symmetrical with respect to w = though more 
spectral weight lies on the electron excitations. S is 
therefore weakly negative over a large temperature range. 
When (f) = 180°, the hole (electron) transmission is high 
(low) [Fig. H(c)]. Since the holes are the main carriers, the 
thermopower is positive. In summary when the quantum 
dot is inserted in the AB interferometer, the magnitude 



and the sign of S can be changed by varying the AB 
phase (/) or magnetic fields threading the AB ring. 

We now address the effect of the energy level 
scheme(ELS) in a quantum dot on the thermopower when 
To = 0. The sign and the magnitude ofthe thermopower 
can be controlled by adjusting the ELS in a dot, too. 
The results in Fig. |^ are computed for the dot with 
^ En+i{U — > 00). In analogy to the heavy 
fcrmion systcms(IIFS), this ELS is equivalent to the Ce 
alloys where the doubly occupied /-electron orbitals lie 
well above the empty /-electron orbital state. Adjust- 
ing the gate voltage capacitatively coupled to the dot, 
the ELS of a dot can be changed in a continuous man- 
ner. When En^i > En+i, more spectral weight of the 
Kondo resonance peak lies on the hole excitations. The 
main carriers are hole excitations, leading to a positive 
thermopower. This inverted ELS is a particle-hole sym- 
metric image of an ELS Ejsi-i < E^^i with respect to 
the point i?7v-i — E^+i and corresponds to the Yb al- 
loys in HFS. Since the flow of electrons is Kondo- assisted 
in a quantum dot, the sign of S is negative (positive) 
when Em-1 < {>)En+i, respectively. In the Ce or Yb 
HFS, electrons are Kondo-scattered so that the sign of S 
is positive (negative) for the Ce(Yb) alloys, respectively. 
Note that the sign of the thermopower is opposite in a 
quantum dot and the bulk HFS with the same energy 
level scheme. When the ELS in a dot is particle-hole 
symmetric(i?Ar_i = E^+i), the Kondo resonance peak of 
the dot is also symmetric with respect to the Fermi level. 
Since the heat currents carried by electrons and holes 
are canceled by each other, the resulting thermopower is 
zero. The linear conductance measures the symmetrical 
part of T(w) near |a;| < fc^T and its value is not sensitive 
to the shape of T(w). On the other hand, the sign of the 
thermopower is sensitive to the degree of the particle- hole 
asymmetry in T(w). Hence the sign of S can probe the 
energy level scheme in a quantum dot. 

In summary, we studied the thermopower in a quantum 
dot in the Kondo regime. In contrast to the bulk heavy 
fcrmion systems, the sign of the thermopower in a quan- 
tum dot can be changed by adjusting the energy level 
scheme or the particle- hole asymmetry in a dot. When 
the dot is inserted in an AB interferometer, the dramatic 
variations in the shape of the transmission spectral func- 
tion T{uj) are possible by controlling the AB phase and 
the tunneling matrices. The rich asymmetrical shapes 
of T{ijj) manifest themselves in the magnitude and the 
sign of the thermopower. Combination of the Kondo ef- 
fect and the Fano interference enables us to control the 
magnitude and the sign of thermopower. 
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FIG. 1: Schematic display of Aharonov-Bohm(AB) interfer- 
ometer with an embedded quantum dot. The magnetic AB 
phase = 27r$ • e/hc is included in the tunneling matrices 
as VdLTLRVRd = \VdLTLRVRd\e"''. * is the magnetic flux 
threading through the AB ring. 
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FIG. 2: Dependence of the transmission probability T{u}) near 
the Fermi level on temperature T, the AB phase (f>. The 
model parameters are chosen as To = 0.5, ed/D = —0.5, and 
Fl = Ffl = 0.07D. With lowering T, the Kondo correlsr 
tion sharpens the shape of T{u}). T is varied as T/Tk = 
100,20,4,0.8,0.16,3.2 x 10"^, 6.3 x 10"^, 1.3 x 10"^. The 
last four temperature curves cannot bo distinguished with the 
naked eye. Fano interference between the direct path and the 
Kondo-resonant tunneling leads to the strong dependence of 
T{ijS) near the Fermi level on the AB phase. 
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FIG. 3: Thermopower S. (i) The solid line is the thermopower 
when the direct tunneling is absent or To — 0. Due to the 
Kondo effect, S is enhanced near the Kondo temperature and 
negative in the energy level scheme En-i i5jv-i-i. (ii) The 
five broken lines show the dependence of S on the AB phase 
(p. Due to Fano interference, the magnitude and sign of 5 are 
modulated by the magnetic AB flux. 



